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o : 

■ Abstract 

o ■ 

■ The effects of noncommutativity in the phase-space of the classical and quantum cosmology of 
Bianchi models are investigated. Exact solutions in both commutative and noncommutative cases 
are presented and compared. Further, the Noether symmetries of the Bianchi class A spacetimes 

\ are studied in both cases and similarities and differences are discussed. 

^ ■ PACS numbers: 04.20.-q, 04.20.Cv 

7— I . 

> ■ 1 Introduction 

m ; 

Since 1947 when noncommutativity between spacetime coordinates was first introduced by Snyder 
t— ( ■ [1], a great deal of interest has been generated in this area of research [2]-[4]. This interest has been 

gathering pace in recent years because of strong motivations in the development of string and M- 
theories, [5, 6]. However, noncommutative theories may also be justified in their own right because 
of the interesting predictions they have made in particle physics, a few example of which are: the 
IR/UV mixing and non- locality [7], Lorentz violation [8] and new physics at very short distance scales 
! Noncommutative versions of ordinary quantum [12] and classical mechanics [13] have also 

been studied and shown to be equivalent to their commutative versions if an external magnetic field 
is added to the Hamiltonian. The behavior of relativistic quantum particles (fermions and bosons) 
^ , in a noncommutative space is investigated in [14]. 

It is a generally accepted practice to introduce noncommutativity either through the coordinates 
or fields, which may be called geometrical or phase-space noncommutativity respectively. Noncom- 
mutative ordinary (quantum) field theories where the geometry is considered as noncommutative 
are interesting to study since they could provide an effective theory bridging the gap between ordi- 
nary quantum field theory and string theory, currently considered as the most important choice for 
quantization of gravity [9]. A different approach to noncommutativity is through the introduction of 
noncommutative fields [15], that is, fields or their conjugate momenta are taken as noncommuting. 
These effective theories can address some of the problems in ordinary quantum field theory e.g. reg- 
ularization [15] and, predict new phenomenon, such as Lorentz violation [16], considered as one of 
the general predictions of quantum gravity theories [17]. 

Since cosmology can test physics at energies that are much higher than those which the experi- 
ments on earth can achieve, it seems natural that the effects of quantum gravity could be observed in 
this context. Therefore, until a completely satisfactory theory regarding cosmology can be afforded 
by string theory, the study of the general properties of quantum gravity through cosmological systems 
such as the universe seems reasonably promising. There are several approaches in considering the 
notion of noncommutativity in cosmology, which, as was mentioned above could be the best alterna- 
tive in the absence of a comprehensive and satisfactory theory from string theory. These ideas have 
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been studied in different works, an example of which can be found in [18, 19]. There, taking coor- 
dinates as noncommuting, it has been shown that noncommutativity affects the spectrum of Cosmic 
Microwave Background [18, 19]. In [18], noncommutative geometry suggests a non-local inflaton field 
that changes the gaussianity and isotropy properties of fluctuations. In this approach the parameter 
of noncommutativity is a dynamical quantity. Introduction of noncommutativity in [19] is based 
on the stringy spacetime uncertainty principle (SSUR), see also [20]. This approach results in the 
appearance of a critical time for each mode at which SSUR is saturated, and which is taken to be the 
time the mode is generated. In this case, in contrast to [18], non-locality appears in time. In both the 
above approaches the spacetime is taken to be the isotropic FRW spacetime. We follow a phase-space 
approach to noncommutativity in cosmology, like that presented in [15] which is different from what 
has been studied in [18, 19] and start with Bianchi models as anisotropic spacetimes rather than an 
isotropic universe. 

In cosmological systems, since the scale factors, matter fields and their conjugate momenta play 
the role of dynamical variables of the system, introduction of noncommutativity by adopting the ap- 
proach discussed above is particularly relevant. The resulting noncommutative classical and quantum 
cosmology of such models have been studied in different works [21]. These and similar works have 
opened a new window through which some of problems related to cosmology can be looked at and, 
hopefully, resolved. For example, an investigation of the cosmological constant problem can be found 
in [22]. In [23] the same problem is carried over to the Kaluza-Klein cosmology. The problem of 
compactification and stabilization of the extra dimensions in multidimensional cosmology may also 
be addressed using noncommutative ideas [24]. 

In this paper we deal with noncommutativity in Bianchi class A models. The classical and quan- 
tum solutions of Bianchi models have been studied by many authors, see for example [25]- [28]. Since 
the Bianchi models have different scale factors in different directions, they are suitable candidates 
for studying noncommutative cosmology. Here, our aim is to introduce noncommutative scale factors 
in Bianchi spacetimes and compare and contrast their solutions to that of the commutative case at 
both the classical and quantum levels. We study the effects of noncommutativity on the underlying 
symmetries of these models [29]-[31] and, show that noncommutativity manifests itself in changing 
the number of such symmetries. It should be emphasized that when we speak of noncommutativity 
in this work, we mean noncommutativity in the fields (scale factors) and not the coordinates, that is 
to say that we study noncommutativity within the context of phase-space only. 



Let us make a quick review of some of the important results in the Bianchi class A models and obtain 
their Lagrangian and Hamiltonian in the ADM decomposition, for more details see [25]. The Bianchi 
models are the most general homogeneous cosmological solutions of the Einstein field equations which 
admit a 3-dimensional isometry group, i.e. their spatially homogeneous sections are invariant under 
the action of a 3-dimensional Lie group. In the Misner [32] notation, the metric of the Bianchi models 
can be written as 



Also, in metric (1), the 1-forms u/ represent the invariant 1-forms of the corresponding isometry 
group and satisfy the following Lie algebra 
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where C l - k are the structure constants. Indeed, the Bianchi models are grouped by their structure 
constants into classes A and B. Because of the difficultly in formulating the class B Bianchi models 
in the context of the ADM decomposition and canonical quantization [33] , it is usually the case that 
one confines attention to the class A models where the structure constants obey the relation Cjj = 0. 
The Einstein-Hilbert action is given by (we work in units where c = K = \6irG = 1) 

S = j d 4 x^d(K - A), (4) 

where g is the determinant of the metric, 1Z is the scalar curvature of the spacetime metric (1) and 
A is the cosmological constant. In terms of the ADM variables, action (4) can be written as [34] 

S = J dtd 3 xC = J dtd 3 xNVh (KijK ij -K 2 + i? - A) , (5) 

where if y are the components of extrinsic curvature (second fundamental form) which represent how 
much the spatial space hij is curved in the way it sits in the spacetime manifold. Also, h and R are 
the determinant and scalar curvature of the spatial geometry hij respectively, and K represents the 
trace of K^j. The extrinsic curvature is given by 

where N^j represents the covariant derivative with respect to hij. Using (1) and (2) we obtain the 
non- vanishing components of the extrinsic curvature and its trace as follows 

(7) 

K 33 = -l(u-2v)e 2 ^~ 2v \ 

* - -4- 

where a dot represents differentiation with respect to t. The scalar curvature R of a spatial hyper- 
surface is a function of v and w and can be written in terms of the structure constants as [27] 

R = C} k C l mn h a h k V + 2C) k Clh? 1 . (8) 

The Lagrangian for the Bianchi class A models may now be written by substituting the above 
results into action (5), giving 



3u 
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+ v 2 + w 2 ) + Ne 3 " (R - A) . (9) 
The momenta conjugate to the dynamical variables are given by 



d£ 12 ■ 3u dC 12 . 3n 8C 12 . 3u 

P - = ^ = ~N Ue Pv = ^=N Ve Pw = dv = N We ' (10) 

leading to the following Hamiltonian 

n = Y A Ne ~* u (- p * +p * + p ™) ~ Ne3u - A ) • ( n ) 

The preliminary set-up for writing the dynamical equations at both the classical and quantum levels 
is now complete. In what follows, we will study these equations in commutative and noncommutative 
cases. 
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3 Classical cosmology 



3.1 Commutative case 

The classical and quantum solutions of the Bianchi models are studied in many works [25]- [28]. Since 
our aim here is to compare these solutions to the solutions of the noncommutative model, in the 
following two sections we consider only the simplest Bianchi class A model, namely type I. The 
structure constants of the Bianchi type I are all zero, that is Cj fc = 0. It then follows from equation 
(8) that R = 0. Thus, with the choice of the harmonic time gauge N = e 3u [35], the Hamiltonian can 
be write as 

H = ±(- P l + pl+pl)+Ae^. (12) 

The Poisson brackets for the classical phase-space variables are 

{xi,Xj} = 0, {pi,Pj} = 0, {xi,p j } = 5 ij , (13) 

where Xi(i = 1,2,3) = u,v,w and pi(i = 1,2,3) = p u ->Pv>Pw Therefore, the equations of motion 
become 

u = {u,H} = -^ Pu , p u = {p u ,H} = -6Ae 6u , (14) 
v = {v,H} = —Pv, Pv = {Pv,T-L} = 0, (15) 

w = {w,H} = y^Pw, Pw = {Pw,7-t} = 0. (16) 
Equations (15) and (16) can be immediately integrated to yield 

Pv = POv, V(t) = J^POvt + V , (17) 

Pw = POw, w(t) = —Powt + w , (18) 

where vo, wq, po v and po w are integrating constants. To integrate equations (14) we note that the 
first integral of this system gives 

Pu = 4A 2 - - A pl (19) 

where A is a constant. Now, integration of this equation depends on the sign of the cosmological 
constant A. If A > 0, then from the second equation of (14) we have p u < 0, and therefore equations 
(19) and (14) result in 

p u {t) =4AcothA(t + t ), (20) 

1 2A 2 r i 
u(t) = -ln—[coth 2 A(t + t )-l\, (21) 

where to is another constant of integration. In the case of a negative cosmological constant, A < 0, 
we have p u > 0, and integration of equations (19) and (14) leads to 

p u (t) = AAtanhA(t + t ), (22) 

1 2A 2 r i 
u{t) = - In [l - tanh 2 A(t + t )\ . (23) 

Finally for a universe with zero cosmological constant, A = 0, we have p u = and the equation of 
motion for u is similar to that of v and w 

Pu = Pou, u(t) = - j^Pout + Uq. (24) 

For a non-zero cosmological constant A these solutions show a singularity at u —> — oo (e 2u — > 0) as 
t — > oo, i.e. the universe evolves to a vanishing spatial size. In the case of a zero cosmological constant 
this behavior depends on the sign of po u . For positive po u we have again a similar singularity as before, 
but for a negative value of po u , the scale factor goes to infinity as t — > oo. While independent of the 
cosmological constant's value, the anisotropic parameters increase to infinity. 
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3.2 Noncommutative case 



Let us now concentrate on the noncommutativity concepts in classical cosmology. Noncommutativity 
in classical physics [13] is described by a deformed product, also known as the Moyal product law 
between two arbitrary functions of position and momenta as 



(/ *q g)(x) = exp 

such that 



f(x 1 )g(x 2 )\ (25) 



aab = I 5ij n Ul3 I, (26) 

\ °ij a ij Pij ) 

where the NxN matrices 6 and f3 are assumed to be antisymmetric with 2N being the dimension of the 
classical phase-space and represent the noncommutativity in coordinates and momenta respectively. 
With this product law, the deformed Poisson brackets can be written as 

{f,g}a = f* a g- g* a f. (27) 

A simple calculation shows that 

%j}a = @ijt \XiiPj\a = $ij \Pii Pj}a = Pij- (28) 

Now, consider the following transformations on the classical phase-space 

x'i = Xi- -Oijpi, p'i =Pi + ^PijX J . (29) 

It can easily be checked that if (xi,pj) obey the usual Poisson algebra (13), then 

{x'i, x'j} = 6ij, {x'^p'j} = Sij + aij, {p'i,p'j} = Pij, (30) 

where = — | {&iPkj + Pi^kj^j- These commutative relations are the same as (28). Consequently, 
for introducing noncommutativity, it is more convenient to work with Poisson brackets (30) than 
a-star deformed Poisson brackets (28). It is important to note that the relations represented by 
equations (28) are defined in the spirit of the Moyal product given above. However, in the relations 
defined by (30), the variables (xi,pj) obey the usual Poisson bracket relations so that the two sets of 
deformed and ordinary Poisson brackets represented by relations (28) and (30) should be considered 
as distinct. 

In this work we consider a noncommutative phase-space in which Pij = and so = 0, i.e. the 
Poisson brackets of the phase space variables are as follows 

{Unc, Vnc} = #1, {v n c, W nc ) = 6 2 , {w nc , U nc ) = 8 3 , 

(31) 

{XinciPjnc} = $ij > {Pinci Pjnc\ = 0. 

With the noncommutative phase-space defined above, we consider the Hamiltonian of the noncom- 
mutative model as having the same functional form as equation (12), but in which the dynamical 
variables satisfy the above deformed Poisson brackets, that is 

U nc = ± (-pl nc + pl nc + p 2 Wnc ) + Ae 6 «- . (32) 
Therefore, the equations of motion read 

U nc = {U nc , H-nc} = -J^Pu nc , Pu nc = {Pu^^nc} = -6Ae 6M " c , (33) 
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Vnc = incline} = yg>v nc ~ QO x ke 6Unc , p Vnc = {p Vnc ,U nc } = 0, 

w nc = {w nc ,H nc } = —p Wnc + 6<9 3 Ae 6 " nc , p Wnc = {p Wnc ,H nc } = 



(34) 
(35) 



Equations (33) are similar to equations (14) in the commutative case. Their solutions are therefore 
as follows 

1 2A 2 



u nc (t) = | In [coth 2 A(t + t ) - 1 



p Unc (t) = 4AcothA(t + t ), 
for positive cosmological constant, and 



(36) 



1 2A 2 r 
u nc (t) = - In [l - tanh 2 A(t + t ) 



P«nc(*) = 4Atanhyl(t + t ), 



(37) 



in the case of a negative cosmological constant. Also for a zero cosmological constant the solutions 
can be written as 



Unc(t) 



1 
12 



POut + UQ, Pu nc =P0u- 



(38) 



Substituting the above results in equations (34) and (35) yields the following equations for v nc (t) and 

W nc {t) 



if A > 0, and 



v nc(t) = —pOvt + 4:6iAcothA(t + to), Pv nc =P0v, 



Wnc{t) = — po w t-46 3 AcothA(t + t ), p Wnc =Po w , 



v nc{t) = j^Povt + 46>iAtanhA(t + t ), p Vnc = p 0v , 



Wnc(t) = -^Powt - 46*3,4 tanh A(t + t ), p Wnc = Pow, 



(39) 



(40) 



if A < 0. Finally if A = 0, the dynamical equations for v(t) and w(t) are again similar to the 
commutative case with solutions 



v nc(t) = J^POvt + VQ, Pv nc = POv, 



Wnc{t) = —POwt + W , p Wnc = P0w 



(41) 



As mentioned before, instead of dealing with the noncommutative variables we can construct, with 
the help of the transformations (29), a set of commutative dynamical variables {u,v,w,p u ,p v ,p w } 
obeying the usual Poisson brackets (13) which, for the problem at hand read 



Pu„.c — Put Pv nc — Pvi Pw nc — Pw 1 

Unc = U- ^9ip v + ^03Pw, 



Vnc = V + ^9ip u ~ ^0 2 p w , 
W nc = W - ^03P U + ^0 2 Pv 



(42) 
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In terms of these commutative variables the Hamiltonian takes the form 

H = ± (- P l+p 2 v +pl)+Ae^-^H°sP™). (43) 
Therefore, we have the following equations of motion 

u = {u,H} = -^ Pu , p u = { Pu ,H} = -QAe^-^+h^), (44) 

v = {v,H} = ^p v -3A0 1 e e ^- 1 2 e ^ +1 2^), Pv = { Pv ,H} = 0, (45) 

w = {w, H} = ^ Pw + ^M^ u -¥^+¥^), Pw = { Pw: H} = 0. (46) 

The solutions of the above equations can be straightforwardly obtained in the same manner as that 
of the system (14)- (16). It is easy to check that the action of transformations (42) on the solutions 
of system (44)-(46) recovers solutions (36)-(41). 

The differences between classical commutative and noncommutative cosmologies are now notable. 
It is clear from solutions in these two cases that both models have a similar singularity at t — > +oo. 
Thus, noncommutativity cannot remove this singular behavior from the commutative solutions. The 
other feature of these solutions is related to their isotropic behavior. With the choice of special 
initial conditions v(0) = p v (0) = and w(0) = p w (0) = for the case A / 0, the commutative 
classical cosmology predicts an isotropic universe, while applying the same initial conditions on the 
noncommutative solutions yields an anisotropic cosmology. 



4 Quantum cosmology 
4.1 Commutative case 

Now, let us quantize the model described above. The quantum cosmology of Bianchi models is well 
studied, see [25]- [28]. Here, for comparison purposes and studying the noncommutativity effects on 
the solutions, we first discuss the commutative quantum cosmology of our model. For this purpose, we 
quantize the dynamical variables of the model with the use of the canonical quantization procedure 
that leads to the Wheeler-DeWitt (WD) equation = 0. Here, H is the operator form of the 
Hamiltonian given by (12) and * is the wave function of the universe, a function of spatial geometry 
and matter fields, if they exist. With the replacement p u — > — and similarly for p v and p w in (12) 
the WD equation reads 



^-^-^ + 24Ae- 
dv? dv 2 dw 2 



i&(u,v,w) = 0. 



(47) 



The solutions of the above differential equation are separable and may be written in the form 
^(u, v,w) = U \u)V (v)W (w) , leading to 



1 d 2 V 
V dv 2 

d 2 U 



1 d 2 W 



du 2 



W dw 2 
+ ( 24Ae 6u ± 9v 2 ) U 



±K 2 , 



o, 



(48) 



(49) 



where r\ and k are separation constants and 9u 2 = rj 2 + k 2 . Equations (48) have simple solutions in 
the form of exponential functions e iziriv (or e ±r)V ) and e =tm "' (or e zizKW ). Also the solutions of equation 
(49) can be written in terms of the Bessel functions as 



U{u) = J v 
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for a positive cosmological constant and 



U(u) = K iv 2 1 



'2|A| su 



for a negative cosmological constant. Thus, the eigenfunctions of the WD equation can be written as 

(50) 



V„(u,v,w) = e~^+™)j„ 2J^e 3u | , A > 0, 




v, w) = e^ v+KW ^K tl/ ( 2 J ^e 3u ] , A < 0, 



(51) 



where for having well defined functions we use the separation constants (48) with plus sign when 
A > and minus sign in the case A < 0. We may now write the general solutions to the WD 
equations as a superposition of the eigenfunctions 

/+oo 
C(u)^i 1/ (u,v,w)du, (52) 

where C(v) can be chosen as a shifted Gaussian weight function e - a ( u - b ) 2 [21]. 
4.2 Noncommutative case 

To study noncommutativity at the quantum level we follow the same procedure as before, namely 
the canonical transition from classical to quantum mechanics by replacing the Poisson brackets with 
the corresponding Dirac commutators as { } — > — i[ ]. Thus, the commutation relations between our 
dynamical variables should be modified as follows 

[UncVnc] = 101, [Vnc,W nc ] = i9 2 , [w nc , U nc ) = l6 3 , [u nc ,p u ] = [v nc ,p v ] = [w nc ,p w ] = %. (53) 

The corresponding WD equation can be obtained by the modification of the operator product in (47) 
with the Moyal deformed product [21] 



-vl + Vl + vl + 24Ae 6u ] * v, w) = 0. 
Using the definition of the Moyal product (25), it may be shown that [21] 

f(u, v, w) * *(u, v, w) = f(u nc , v nc , w nc )^(u, v, w), 



(54) 



(55) 



where the relations between the noncommutative variables u nc , v nc , w nc and commutative variables 
u, v, w are given by (42). Therefore, the noncommutative version of the WD equation can be written 

as 

- Q2 Q2 Q2 



_ _ _T_ + 24Ae 6(M ~5 e iP-+l e 3P-) 

dv? dv 2 dw 2 



W(u,v,w) = 0. 

We again separate the solutions into the form ^(u,v,w) = e l ^ v+KW ^U(u). Noting that 



(56) 



e 6(«-^ij*+^ 3 p»)gr( UjUjU; ) = e ^( UjV + 3i0i,w-3i6 3 ) 

— e 6ti^j^ u ^ e i7/(v+3z0i) e iK(tu-3ie3) 

_ e 6«f/^ n ^ e «(»?f+Kw) e -3(»;ei-K6»3) 

= e 6u e- 3(Tiei ^ Kds) ^(u,v,w), 



(57) 
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Figure 1: The figure on the left shows the square of the commutative wavefunction while the figure on the right, the 
square of the noncommutative wavefunction. We choose w = for simplicity. 



equation (56) takes the form 

- Q2 Q2 Q2 



du 2 dv 2 dw 2 



V(u,v,w) =0. (58) 



Following the solutions to equation (47), we can find the general solutions for the above equation as 

r+oo 

y(u,v,w) = / A(u)^ v (u,v,w)du, (59) 



where again, A{y) is chosen as a shifted Gaussian function and ^f u (u, v,w) are the eigenfunctions of 
equation (58), that is 

v, w) = e~^ v+KW \j u ( 2W ^ e 3(«-|^ 1+ i K e 3 ) ] > A > 0, (60) 



u, w) = e i{riv+Kw) K iu ( 2y ^ e 3 («-|^i+l K0 3)j ? A < 0. (61) 

where again 9z^ 2 = rj 2 + k 2 . Figure 1 shows the square of wave functions of the commutative and 
noncommutative universes. The comparison of the figures show that the noncommutativity causes 
a shift in the minimum of the values of u, corresponding to the spatial volume. The emergence of 
new peaks in the noncommutative wave packet may be interpreted as a representation of different 
quantum states that may communicate with each other through tunnelling. This means that there are 
different possible universes (states) from which our present universe could have evolved and tunnelled 
in the past, from one universe (state) to another (see the first reference in [21]). 

5 Symmetries of commutative and noncommutative Bianchi class 
A spacetimes 

In this section we investigate the Noether Symmetries of the Bianchi class A models and study the 
effects of noncommutativity on these symmetries. Following [29] and [30], we define the Noether sym- 
metry of the spacetime as a vector field X on the tangent space of the phase-space-like (u, v, w, u, v, w) 
through 

X ® + Q ® + d ^ da d ^ d(3 d ^ d 
a du dv ^ dw dt dii dt dv dt dw ' 
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such that the Lie derivative of the Lagrangian with respect to this vector field vanishes 

L X C = 0. (63) 

In (62), (a,/3,j) are linear functions 1 of (u,v,w) and ^ represents the Lie derivative along the 
dynamical vector field, that is 

d . d . d . d 



dt U du V dv W dw 



(64) 



It is easy to find the constants of motion corresponding to such symmetry. Indeed, equation (63) can 
be rewritten as 



' a du dt dii J V dv dt dv J \ dw dt dw , 



(65) 



Noting that §| = we have 



which yields 



dp u da \ ( dp v d(3 \ ( dp w d-y 



d 



o. 



— {ap u + (3p v + -fp w ) = 0. 



Thus the constants of motion are found as 

Q = ap u + (3p v + 7P U 



(66) 
(67) 
(68) 



In order to obtain the functions a, (3 and 7 we use equation (65). However, since in the noncommuta- 
tive case we use the Hamiltonian formalism, equation (66) is better suited in finding these coefficients 
which, equivalently, can be written as 



a{Pu, H} + P{ Pv ,H} + 7{Pw,n} + 

Pv 



da , „ da , „ da , „ 



9I} , w^ 8 ! 3 , ux^Ws ux 



|2 { „. w} + |2 { „, w} + |l {liI , w} 



Pu + 

p w = 0. (69) 



In general, the expression above gives a quadratic polynomial in terms of momenta with coefficients 
being partial derivatives of a, [3 and 7 with respect to the configuration variables u, v and w. Thus, 
the expression is identically equal to zero if and only if these coefficients are zero, leading to a system 
of partial differential equations for a, j3 and 7. In the following subsections we obtain such symmetries 
for Bianchi class A models in the commutative and noncommutative cases. 



5.1 The Bianchi type I model 

This model is the simplest in that all the structure constants are zero. The Hamiltonian in the gauge 
N = e 3u is given by equation (12). Substituting the Poisson brackets (14)- (16) in equation (69), we 
get 

. 1 da n 1 d(3 o 1 c?7 2 I f da d(3\ 1 f da dj\ 

- SA " e - 12 + 12 + l28^» + 12 { Si - Si) + 12 [to - Yu) P " P " 

^(S + S)^^ ' <7o) 



1 One may, of course, seek more complicated symmetries. However, we have chosen the linear case for simplicity. 
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which leads to the following system 



da d(3 d^/ 

a = °' 7r = 7T = 7T 1 = ' 
ou ov dw 



da 8(3 da d~/ d[3 dj 
dv du dw du dw dv 



(71) 



The above system has the general solutions [29, 30] 

a = 0, (3 = aw + b, 7 = -av + b' , (72) 
where a, b and b' are constants. Thus, according to (68) the independent constants of motion are 

Qi = Pv, Q2 = Pw, Q3 = wPv ~ vp w , (73) 
which correspond to the following symmetries 

d d d d d d 

Xi = — , X 2 = —, X 3 =w 7) v-= Vw-^r -v-z-r. (74) 

dv dw dv dw dv dw 

These symmetries satisfy the Lie algebra 

[X 1 ,X 2 ]=0, [X 2 ,X 3 ]=X 1 , [X 3 ,X 1 ]=X 2 . (75) 

Now, let us see which of the above symmetries survive in the noncommutative case. Here, the 
Hamiltonian is given by (43) and the required Poisson brackets are given by equations (44)- (46). 
Substituting these Poisson brackets into relation (69) gives 

12 du Pu + 12 dv Pv + 12 dw Pw + 
1 (da dfi\ 1 (da 9 7 \ 1 (8(3 d-/\ 

uldv-dn PuPv + 12 ( ^ - m J PuPw + 12 + dv j PvPw + 



, da n da\ / d(3 n dfi\ / dj n &y\ 
dv dw J V dv dw / V ov dw . 



0. (76) 



Thus, in addition to equations (71) we have a number of extra restrictions on a, (3 and 7 as 

, 1 ^- h ^ = ei f-,^ = 91 §2- % |I = 0, (77) 
dv dw dv dw dv dw 

which should hold for all values of 9\ and #3. Imposing relations (77) on the solutions (71) results in 
(3 and 7 becoming constants. Therefore, in the noncommutative case the only constants of motions 
are 

(Ql)nc=Pv, (Q2)nc=Pw, (78) 

which correspond to the symmetries 

(X 1 ) nc =-^, (X 2 ) nc = —, [X 1} X 2 ]=0. (79) 

Consequently, the third symmetry is absent in the noncommutative case. A quick look at the Hamil- 
tonians (12) and (43) shows that the variables v and w are cyclic and consequently the corresponding 
momenta are constants of motion. In addition, Hamiltonian (12) is rotationally invariant about the 
u-axis, that is v 2 + w 2 . Thus, the angular momentum about this axis is conserved. This symmetry 
is absent in the noncommutative Hamiltonian (43). 

At this point, it is appropriate to investigate the symmetries when the cosmological constant is 
zero. In this case there is no difference between commutative and noncommutative Bianchi type I 
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models. It is clear from (70) that the symmetries in this case are obtained from the system (71) 
without the constraint a = and that the solutions are given by 



a = av 



+ bw + c, (3 = au + b'w + c , 7 = bu — b'v + c" , (80) 



with a, b, V , c, d and c" being constants. After a simple calculation we find the independent constants 
of motion, given by 



(81) 



(82) 



Ql=Pu, Q2=Pv, Qz=Pw, 

Qi = wp v -vp w , Q 5 = wp u + up w , Q 6 = vp u + up v , 
which correspond to the following symmetries 

d d d 

X\ = — , X 2 — — , X 3 — — , 
ou av aw 

d d d d 
X 4 = w- v— + w— - v— 

av aw av aw 

d d . d . d 

ou aw ou aw 

d d . d . d 
Xr = v— — h u— — \- v— — h«— — , 

ait aw au av 

satisfying the following Lie algebra 

LY 1? X 2 ] = 0, [X 1 ,X 3 ] = 0, [X 1; X 4 ] = 0, 

[Ai,X 5 ] = X 3 , [Xi,X 6 ]=X 2 , [X 2 ,X 3 ]=0, 

[X 2 ,X 4 ] = -X 3 , [X 2 ,X 5 ]=0, [X 2 ,X 6 ] = X 1 , (83) 

[^3,-^4] = X 2 , [X 3 ,X 5 ]=Xi, [X s ,Xq]=0, 

[X a ,X 5 ] = -Xq, [X 4 ,X 6 ]=X 5 , [X 5 ,Xq] = X 4i . 

5.2 The Bianchi type II model 

This model is characterized by the 1-forms 

uo 1 = dx — zdy, uj 1 = dy, u> 3 = dz. (84) 
Therefore, the non-vanishing structure constants may be obtained from equation (3) 

C23 = _C 32 = 1> ( 85 ) 
which, upon substitution into relation (8), yields the scalar curvature of the corresponding 3-geometry 

R = 2e ~ 2u+Av+A ^ w . (86) 

Thus, the Hamiltonian for the model can be written from equation (11) with the result (again we use 
the gauge N = e 3u ) 

H = -L (-pi + pi + p2 ) + A e 6« _ 2e 4 («+-+v / 3») . ( 87 ) 

In order to obtain the symmetries, we have to solve equation (69) for a, (3 and 7. The required 
Poisson brackets are 



as 



{Pu,T-C} = - 



6Ae 6n - 8e 4{u+v+V3w) 



{ Pv ,H} = 8e 4 ( u+v+V3 ~ w \ { Pw ,H} = 8V3e^ u+v+ ^ w \ 
{u,H} = -j^Pu, {v,H} = — p v , {w,H} = —Pw- 
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These Poisson brackets make equation (69) read 

+ 5 i Q + p + V^7je 12du Pu 12dv Pv 12dw Pw 

T2 - I) ^ + L2 " ID ^ + I2 + p ^ = °' (89) 

which, in turn, results in the following system 



i- da 80 8j 

a = 0, a + /3 + V3 7 = 0, — = -£ = -1=0, 

era 8v ow 



da 80 q da d^f ^ 80 d^y ^ 
8v 8u 8w 8u 8w 8v 



(90) 



The solutions of this system are similar to the solutions of the system (71) with the additional 
constraint a + + \/3~7 = 0. Therefore, we obtain a = 0, and 7 = const, as solutions, i.e. we have 
a = 0, = — y/3 and 7 = 1. Consequently, there is only one symmetry 

^ 8v <9u>' ^ 

with the corresponding constant of motion 

Q = -V3 Pv +p w . (92) 

In the noncommutative case we have the system of equations (90) plus the additional equations in 
terms of the noncommutativity parameters and partial derivative of a, and 7, similar to the last 
term in (76). It is easy to see that these additional equations do not have any effect on the solutions 
of system (90), and thus, also hold in the noncommutative phase-space. 

Now, let us consider this model with a zero cosmological constant. In this case the restriction 
a = is removed from the system of equations (90) and it may be shown that the general solutions 
can be written as 

a = VSav — aw — c — V3c, = VSau + aw + c, 7 = — au — av + c , (93) 
where a, c and d are constant. Therefore, in this case we have three independent constants of motion 

Qi = -V3p u + p w , Q 2 = -p u + Pv, 

(94) 

Q3 = (VSv - w^j p u + (VSu + w^j p v - (u + v)p w , 

which correspond to the following symmetries 

/rrd 8 8 8 

Ai = -V3— + — , X 2 = -— + —, 
8u 8w 8u 8v 

X 3 = (VSv - w) £ + (VSu + w) ^ - („ + t,) A + (95) 

satisfying the Lie algebra given by 

[X 1 ,X 2 ]=0, [X 1 ,X 3 ]=V3X 1 -2X 2 , [X 2 ,X 3 ] = -V3X 2 . (96) 
The noncommutative phase-space with zero cosmological constant has the Hamiltonian 

Jl = }_ ^_p2 +p 2 + p^j _ 2e 4 «-2(?ip t ,+26»3P» e 4? ) +26»ip u -26»2P» e v / 3(4« I -26»3p It +202P t ,) i (97) 
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To construct equation (69) we need the following Poisson brackets 

{Pu,H} = { Pv ,H} = ^{ Pw ,H} = 8e A e B e c , 



{u,H} = --^pu - 4 (0i - V36 3 ) e A e B e 
{v,H} = -^Pv + 4 (<?i - V39 3 ) e A e B e 



12'"" ' ' 

(98) 

12 iJ " ^ ' " *° 

{«;, = ^ Pw - 4 (0 3 - 2 ) e A e B e c , 

where A = 4u-29 1 p v + 20 3 p w , B = 4v + 29 1 p u -29 2 p w and C = \/3(4w-29 3 p u + 29 2 p v ). Substitution 
of these relations into equation (69) shows that solutions (93) should satisfy the following additional 
constraints 

( ft _V3fc)§![ + (*-fc)|| = 0, (99) 

^7 _ ^7 = 
<9u dv 

Imposing these conditions on solutions (93), we are led to 

a (V30i - 2 - 20 3 ) = 0. (100) 

Now if a = 0, i.e. a, (3 and 7 are constant, then solutions (93) automatically satisfy conditions (99). 
Therefore, the noncommutative case has the symmetries X\ and X 2 only, given by (95). On the other 
hand if a / 0, there is a special combination of 0's, that is 

v/30i - 2 - 20 3 = 0, (101) 

for which all the symmetries (95) in the commutative case can also be recovered in the noncommu- 
tative case. This phenomena can be understood by noting that in this case Hamiltonian (97) takes 
the form 

H=± (_ p 2 +p 2 + p 2\ _ 2e 2(e 1 -Vse 3 )(2Q 2 -V3Q l ) e 4(u + v+V3w)_ (1Q2) 

Since Q\ and Q 2 are constants of motion the above Hamiltonian differs from the Hamiltonian (87) 
(for A = 0) only by a constant factor in the potential term. As it is clear from equation (69) which 
is the starting point in our procedure for finding the symmetries, this constant factor does not have 
any effect on the symmetries. Thus in the special case (101), Hamiltonians (102) and (87) have the 
same structure as far as the symmetry considerations are concerned. 

5.3 The Bianchi types VI and VII 

These types of Bianchi models are characterized by the 1-forms 

uj 1 = cosh zdx =F sinh zdy, J 2 = — sinh zdx + cosh zdy, uj 3 = dz, (103) 

with the corresponding structure constants 

= ~C\ 2 = ±1, Cl = -C 2 13 = -1, (104) 

where the upper and lower signs denotes the Bianchi type VIo and VIIo respectively. The scalar 
curvature of the corresponding 3-geometry can be evaluated from equation (8) with the result 

R = 4e u+4v [cosh(4v / 3w) ± l] , (105) 
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leading to the Hamiltonian 

To find the symmetries, we need the following Poisson brackets, obtained from relation (69) 

{ Pui H} = -6Ae 6 " + \Qe A{u+v) [cosh(4V3w) ± l] , 
{p v ,H} = 16e 4( - u+v) [cosh(4V3u;) ± l] , 

{ Pw ,H} = 16V3e 4( - u+v hmh{4:V3w), 
{u,H} = -j^Pu, {v,H} = — p v , {w,H} = —p w . 

Substitution of the above Poisson brackets into equation (69) leads to the following equation 



(106) 



(107) 



-6aAe 6u + 16(a + (3)e 4{u+v) [cosh(4v / 3w) ± 1 

1 dp __ 2 1 57 2 



+ 



16V3-fe^ u+v ^ smhUV3w) - — ^pl + 

12 du 



l2 8v Pv+ \2dw Pw + 



1 /da dp\ 
12 \~dv~ ~ du) 



1 (da dj\ 1 (d(3 dj 

PuPv + T77 o TT - ] PuPw + T7: T7- + T7- PdP 



12 \<9ti; <9u 



12 V<9if 9f 



0, 



(108) 



which for a non-zero cosmological constant immediately yields a = P = 7 = 0, i.e. Bianchi types 
VIo and VIIo universes with non-zero cosmological constants have no symmetries and, therefore, 
no constants of motion. It is clear that this also holds when we consider such spacetimes in a 
noncommutative scenario. 

Let us now assume that the cosmological constant in these models is equal to zero. In this case, 
from equation (108) we have 



da dp d^i 
7 = 0, a + P = 0, — = -f = -L = 0, 

du dv dw 



(109) 



da dp _ q da ^7 _ g &P _j_ ^7 
dv du dw du dw dv 



0, 



which admit the solution a = — P = const. Therefore, the only symmetry and its corresponding 
constant of motion are obtained as 

d_ _ d_ 

du dv ' 

The Hamiltonian in the noncommutative model for A = is 



X 



Q =Pu ~Pv 



(110) 



_ Ae {4u-28 1Pv +2e 3 p w ) e (±v+2e iP u-2e 2 p w ) \ cosh ^^ w _ 2V39 3Pu + 2V39 2 p v ) ± 1 
yielding the following Poisson brackets 

{ Pu ,H} = {p v ,H} = 16e A+B (coshCi 1) , {p w ,H} = lev^sinhC, 



(111) 



{u,H} = 

{v,H} = 
{w,H} = 



~T2 Pu ~ 89ieA+B ( coshC± !) + 8 #3 V3e A+B sinhC, 



(112) 



12 
1 

\2 



p v + 86 ie A+B (coshC ± 1) - 86 2 V3e A+B sinhC, 



p w ~ 8(6 3 - 9 2 )e A+B (cosh C ± 1) . 
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Using these relations in equation (69), we conclude that the solutions of the system (109) should also 
satisfy the following additional conditions 



da 


da 


= o, 


a da 




dv 


du 


du 


dv 


dp 


dp 


= o, 


d d_p__ 




dv 


du 


du 


dv 



o, W.-«£-o, 



0, («3-«£=0. 



(113) 



If 02 ^ these equations have the previous solutions a = — P = const., and thus the spacetime has 
the same symmetry as given by (110). However, in the case where 02 = 63, the above system admits 
the solutions 

a = -P = aw + b, (114) 
with constants a and b. The independent constants of motions are then given by 

Qi=Pu~Pv, Q2 = w(.Pu- Pv), (115) 
corresponding to the following symmetries 

d d „ ( d d\ . ( d d 



• Y '=*-*- • Y2 = , °(*-w + n»-W' [Xi - Y2l=a <ll6) 

Here, we have an additional constant of motion relative to the commutative case. This is notable in 
Hamiltonian (111) in that when 62 = 03, it takes the form 

H = ij- (-pi +p 2 v + pi) - 4e 2eiQl e 4 ^ [cosh(4v^ - 20 2 VSQi) ± l] . (117) 

The above Hamiltonian differs from Hamiltonian (106) (for A = 0) not only by a constant factor in 
the potential term (which has no effect on the number of symmetries like the Bianchi type II case) 
but also in the argument of the "cosh" term. This argument is modified by a constant shift which 
results in a new symmetry (compare (115) with (110)). 

5.4 The Bianchi types VIII and IX 

In the Bianchi type VIII model the invariant 1-forms are 

uj 1 = cosh y cos zdx — sin zdy, 

J 2 = coshy sin zdx + cos zdy, (118) 
a; 3 = sinhydx + dz, 

and thus the corresponding non-vanishing structure constants can be evaluated from (3) with the 
result 

C23 = —C32 = —1, C 31 = — Cf 3 = —1, Cf 2 = —C21 = 1- (119) 

The scalar curvature of the spatial hypersurface can be obtained by substituting the above structure 
constants in (8), yielding 

R = -4 e -2u+iv cogh^^) _ 2e -2u~8v + Ae -2u-2v cosh ^V3w) - 2 e - 2u+4v . (120) 
Therefore, we can write the Hamiltonian using (12) 

H = ^ (-pi + p 2 v + p 2 w ) +ke %u + 4e i{u+v ^ cosh(4 v / 3u») - Ae 2 ^^ cosh(2 v / 3u;) + 2e 4{u - 2 '^ +2e< u+v l 

(121) 
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As it is clear from the above form, the spacetime has no symmetry. Indeed, evaluating the required 
Poisson brackets and constructing equation (69), we see that this system admits only the trivial 
solution a = (3 = 7 = 0, i.e. there are no symmetries neither in the commutative nor noncommutative 
models. 

Our final goal is to address the same issues in the Bianchi type IX model, characterized by the 
structure constants 

c jk = £ ijk- (122) 
Calculations similar to the previous models lead to the following Hamiltonian 

U = JL (-pi + p2 + p ty + Ae 6u + e 4(u - 2v) - Ae 2{2u ~ v) cosh(2 v / 3^) + 2e 4{u+v) [cosh(4v / 3w) - 1 . 

(123) 

Again, a glance at the equations resulting from the above Hamiltonian shows that this spacetime has 
no symmetry, neither in the commutative nor in its noncommutative version, even if the cosmological 
constant is equal to zero. 



6 Conclusions 

In this paper we have studied the Bianchi type I classical and quantum cosmology in both commu- 
tative and noncommutative scenarios. We have obtained exact solutions of the vacuum gravitational 
field equations in three cases when the cosmological constant is positive, negative or zero. The cor- 
responding classical cosmology shows a singularity at t — ► 00, together with increasing anisotropic 
parameters v and w. We have seen that noncommutativity can not exclude such behavior from the 
solutions. The corresponding quantum cosmology of the model in both commutative and noncom- 
mutative cases is obtained from the exact solutions of the WD equation, showing how such solutions 
differ. 

We have also studied the Noether symmetries of the Bianchi class A models and investigated their 
noncommutative counterparts. In Bianchi type I, we have shown that for a non-zero cosmological 
constant, there are three Noether symmetries, two of which retain their character in the noncom- 
mutative case. For a zero cosmological constant, there is no difference between commutative and 
noncommutative cases and the number of symmetries is six. The constants of motion in this case are 
the momenta in u, v and w directions and their corresponding angular momenta. However, we note 
that because of the wrong sign (minus) of p u in the Hamiltonians, this dynamical variable behaves 
like a ghost field and thus the angular momenta-like constants Q5 and Qq do not have the usual sign. 
The Bianchi type II model with a non-zero cosmological constant in the commutative and noncom- 
mutative cases has only one symmetry, which is a Killing vector field. If the cosmological constant 
is zero, this spacetime has three symmetries when studied in its noncommutative form. We have 
shown that when the phase-space has a noncommutative structure, one of its symmetries is removed 
except when the noncommutative parameters satisfy a certain relation. With this relation between 
the noncommutative parameters, although the Hamiltonians are different in the commutative and 
noncommutative cases, they show the same structure and thus exhibit the same symmetries. 

The Bianchi types VIo and VIIo with non-zero cosmological constants have no symmetries neither 
in the commutative nor noncommutative versions. But in the case of a zero cosmological constant, 
these spacetimes exhibit one symmetry if their phase-space is commutative. This symmetry remains 
unchanged in the noncommutative phase-space as well. It is interesting to note that if the noncom- 
mutative parameters 62 and #3 are equal, these models have an additional symmetry. In this case 
the corresponding Hamiltonians of the commutative and noncommutative phase-space have the same 
structure but differ in a shift term in the variable w. In the last section of the paper we have dealt 
with the Bianchi type VIII and IX models and shown that they have no symmetries at all. 
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